Phase coexistence in consolidating porous media 
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The appearence of the fluid-rich phase in saturated porous media under the efi'ect of an external 
pressure is investigated. For this purpose we introduce a two field second gradient model allowing 
the complete description of the phenomenon. We study the coexistence profile between poor and 
rich fiuid phases and we show that for a suitable choice of the parameters non-monotonic interfaces 
show up at coexistence. 
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I. INTRODUCTION 

When a porous medium is plunged in an infinite fluid 
reservoir, the solid matrix absorbs fluid until an equi- 
librium state is reached. Many interesting features of 
this swelling phenomenon have been demonstrated exper- 
imentally. The amount of swelled fluid can be controlled 
via different external parameters such as the fluid pres- 
sure in the reservoiri^, that is via its chemical potential, 
the fluid velocity^, or a mechanical pressure exerted on 
the solidi-^. 

The solid-fluid segregation in consolidation is seen 
when, depending on the external pressure acting on the 
porous material, phases differing in fluid content are ob- 
served. This problem has been addressed^ by the authors 
in the framework of continuum mechanics adopting in 
particular a first gradient model; the existence of two 
different phases depending on the external pressure has 
been proved. In that context the interesting question of 
the coexistence of the two phases could not be posed for 
the first gradient nature of the model. In this paper we 
propose a more general two field one dimensional second 
gradient model to study the profiles connecting two co- 
existing phases and the formation of critical droplets, if 
any, of one phase into the other. 

The adopted approach in modeling the behavior of 
porous continua is essentially based on a pure solid La- 
grangian description of motion, referring kinematics to 
the reference configuration of the porous skeleton (see 
Section [ni- 

The constitutive model is purely phenomenological, 
which means that the overall potential energy, regarded 
as a function of the strain of the skeleton and the fluid 
mass density (per unit volume, in the solid reference con- 
figuration), is built up in such a way to describe the ex- 
istence of two states of equilibrium: the solid-rich and 
the fiuid-rich phase. Thus no refined description of solid 
grain connectivity, as well as connection among regions 
with different porosity is available in such a model. Con- 
versely the constitutive state parameters are selected so 
as to describe the showing up of the fluid-rich phase, 
which is possibly associated to the occurrence of fluid 
segregation. 



II. THE POROMECHANICS SETUP 

Let C M be the reference configurations for the solid 
and fluid components^. The solid placement is a C^- 
diffeomorphism Xs{-,t) ■ Bg M. such that Xsi^sjt) 
is the position occupied at time t by the solid particle 
Xs in the reference configuration B^. Consider— (/>{■, t) : 
Bg ^ M. such that (f>{Xs,t) is the fluid particle which 
at time t occupies the same position of the solid parti- 
cle Xg. Assume also 4>{-,t) to be a C^-diffeomorphism, 
thus the map (t>{-,t) associate univocally a solid particle 
to a fluid one and vice versa. The fluid placement map 
Xf(-, t) : M K, giving the position of a fluid particle Xf, 
is defined as xii^iit) •= Xs{'P~^{Xf,t)-it). The current 
configuration Xs{Bs,t) at time t is the set of positions of 
the superposed solid and fluid particles. 

Let Js{Xs,t) := \dxs{Xs,t)/dXs\ be the Jacobian of 
the placement map Xsi',t) measuring the ratio between 
current and reference volumes of the solid component; 
we let e{Xs,t) := (Js(Xs,t)^ — l)/2 be the strain field. 
Let go^a{Xa) with a ~ s, f, be the solid and fluid ref- 
erence densities; we define the fluid mass density field 
■m{Xs,t) := goj{(j>{Xs,t))d(j){Xs,t)/dXs. Assuming that 
the mass is conserved, it is not difficult to proved that 
the field m can be interpreted as the fluid mass density 
measured w.r.t. (with respect to) the solid reference vol- 
ume. 

Assume, now, that the Lagrangian density 
A" 



form 



, Xs,<^) of the system is in the 



if = T(xs>,Xs 



1>(x">",Xs>',Xs,0) (1) 



where T is the kinetic energy density and <& is the overall 
potential energy density accounting for both the internal 
and the external forces. In ([T]) we have denoted with the 
dot the derivative taken w.r.t. time and with the prime 
the derivative w.r.t. the solid reference space variable. 
The equation of motion for the two fields Xs and (j) can be 
derived assuming that the possible motions of the system 
in an interval of time (^1,^2) C M are those such that the 
fields Xs and (j) are extremals for the action functional 



dXg 



dt^(Xs,...,0) (2) 
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in correspondence of the independent variations of the 
two fields Xs and (p on Bs x (^1,^2)- In other words any 
possible motion of the system in the considered interval 
is a solution of the Euler-Lagrange equations associated 
to the variational principle SA = 0. 

If one is interested to find equilibrium profiles XsiXg) 
and 4){Xs) of the system, namely, the solutions of the 
equations of motion independent of time, since the kinetic 
energy associated to those profiles is equal to zero, the 
Lagrangian density reduces to minus the potential energy 
density. In this case the action is given by ^ where the 
time integral gives a not essential multiplicative constant, 
and the variational principle associated to such an action 
gives the seeked for equilibrium profiles. 

Assume, now, that the effect of the internal forces ex- 
changed by the solid and fiuid particles and that of the 
conservative external fields can be described via a po- 
tential energy density $(m',£',m,e) depending on the 
kinematic fields Xs and </> only through the strain and 
the fluid mass density fields. Note that the strain de- 
pends only on Xs and the fluid mass density only on 0, 
hence the independent variations of those primitive fields 
reflect on independent variations of e and m. Thus, lim- 
iting the study to boundary value problems expressed in 
terms of the fields e and m, we can treat the fields e and 
m as primitive, consider their independent variations and 
look for the equilibrium profiles £{Xs) and m{Xs) start- 
ing from the variational principle 

6 f dX, $(m'(X,), £'(X,), m{X,), e{Xs)) = (3) 



Finally, we can derive, starting from ([3]), the equations 
governing the equilibrium profiles e and m. By comput- 
ing the variation of the action functional on Bg = (£1, €2), 
with £i,£2 G we get the Euler-Lagrange equations 



a$ d a$ 

'de ^ dXs de' 



and 



dm 



d a$ 



dXs dm' 



with boundary conditions ensuring that 



+ —5m 
de' dm' 



= 



(4) 



(5) 



where de and 5m are, respectively, the variations of 
the strain and fluid mass density fields. For instance 
Dirichelet boundary conditions would do the job, since we 
would have fe(^i) = 5e{t2) = and 5m{ii) = 5m{l2) = 
0. But for potential energy densities $ at least quadratic 
in the derivatives e' and m' even Neumann boundary 
conditions would be acceptable. In the sequel we shall 
refer to any solution of the Euler-Lagrange equations (|4]) 
with suitable boundary conditions as an equilibrium pro- 
file of the system corresponding to the chosen boundary 
conditions. 

The model ([3|) is called a first gradient model if the 
potential energy density $ does not depend on the first 
derivatives of the strain and of the liquid density, oth- 
erwise the model is said a second gradient model. This 



way of classifying the models is related to the fact that 
both £ and m depends, by definition, on the gradient of 
the primitive kinetic fields Xs and </>. 

Second gradient theories are suitable to be devel- 
oped for modeling stress/strain concentration due, for in- 
stance, to the presence of geometrical singularities (crack 
propagation in fracture mechanics^) or phase transitions 
as in the case of wettingiSi^i. In particular second gradi- 
ent poromechanics has been recently formulated^ii^ ex- 
tending the standard arguments of the Biot theorjsi^. 
Such a model addresses the description of those defor- 
mation phenomena which occur at the same length scale 
as that where high gradients in deformation can be de- 
tected. Classical poromechanicai is not able to describe 
these phenomena: the macroscopic model is regarded 
in that case as the average of a microscopic one where 
a kind of stationarity assumptions^ (spatial crgodicity) 
on the random field which characterizes the microscopic 
mechanical properties of the material has been formu- 
lated. This allows for replacing ensemble averages with 
volume averages insofar as the characteristic size of the 
heterogeneities is much smaller than the typical length 
scale of the reference volume element (RVE). If this is 
no more the case, the classical assumptions of uniform 
strain (stress) or periodic boundary conditions, for ev- 
ery reference volume, are no more valid, but, conversely, 
macroscopic strain gradient plays a crucial role in speci- 
fying the state of stress/strain inside the RVE itself. 

The goal, here, is to formulate a second gradient 
poromcchanical model for describing the transition from 
the standard Biot-like equilibrium, associated to a com- 
pacted solid, versus the fluid-segregated phase describing 
duct thinnering in the matrix and, consequently, fluid 
mass concentration in the pores. Thus, we consider a 
model with total potential energy density in the form 



<I>(m', e' , m, e) = K{m' , e' , m, e) + ^'(to, £) 



(6) 



where K is a polynomial quadratic function of e' and m', 
'I' is a differentiable function positively diverging along 
any radial direction in the plane e-m, having at least a 
local minimum, and whose stationary points arc isolated. 
Since K is quadratic, we have that a constant solution of 
the Euler-Lagrange problem (|4]) must necessarily satisfy 
the equations 'I'^ = and ^'m = 0; in other words a 
constant proflle must be constantly equal to an extremal 
point of the first gradient part 5* of the total potential 
energy. 

We then let a phase of the model to be a constant 
equilibrium profile equal to one of the local minima of 
the function vj/. 

Note that the Euler-Lagrange problem (|4]) and (O for 
the first gradient model associated to namely the one 
obtained for K = 0, is the system of algebraic equations 
'i'm = and ^'^ = 0. Since the stationary points of 
the two variable function are isolated, we have that 
the equilibrium profiles for such a model are necessarily 
constant functions of Xg € Bg equal to the values of the 
stationary points of 
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Hence, in the case of a first gradient model it is not 
possible to discuss phase coexistence, since there exist 
only continuous constant equilibrium profiles. On the 
other hand, in second gradient models, different (not 
constant) equilibrium profiles can exist. This fact al- 
lows us to pose the problem of the coexistence of two 
existing phases. Suppose that the model exhibits the 
two phases (rni,ei) and (m2,£2); a connection^ hctwccn 
those phases is an equilibrium profile m, e of the action 
functional on Bg = (—00, +00) satisfying the boundary 
conditions to(— 00) = ?7ii, m(+oo) = 7712, £(—00) = £1, 
and £(+00) = £2. We say that the two considered phases 
coexist if and only if a connection does exist. 



e, m 
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III. THE MODEL 



FIG. 1: From the bottom to the top the graphs of e^{p), 
ms(p) = bes{p), and mt(p) for a = 0.5, 6=1, and a = 100. 



We study, now, a particular poroelastic model and in that 
framework we discuss the existence of the consolidation 
phase transition and prove the coexistence of the two 
phases for a particular value of the external pressure. 
More precisely we consider the poroelastic system with 
overall potential energy density ([6]) with 

K{m',e') := ^[ki{e'f + 2k2e'm' + fcsK)'] (7) 
with ki, ^3 > 0, ^2 € M such that fcifcs — fc| > 0, and 
^'(m, e,p) := ^m^(3m^-8&£m+66^£^)+*B(m, e,p) (8) 



where 



(9) 



is the Biot potential energy density^, a > is the ratio 
between the fluid and the solid rigidity, 6 > is a cou- 
pling between the fluid and the solid component, p > is 
the external pressure, and a > is a material parameter 
responsible for the showing up of the additional equilib- 
rium. We remark that the condition fcifca — fc| > ensures 
that the second gradient part K of the overall potential 
energy density is convex. Under this assumption there 
exists a minimizer for the action functional 



dXs$(m',£',m,£) 



on a bounded domain. As we will see later to ensure the 
existence of a connection profile, which is a Dirichelet 
problem on an unbounded domain, it will be necessary 
to assume fcifca — fc| > 0; in the limiting case fcifcs — fc| ~ 
the existence of the connection will depend on the choice 
of the parameter fci, ^2, and ^3. 

We have already studied^ the associated first gradient 
model with overall potential energy density 'J' and we 
have proven the existence of a phase transition driven 
by the external pressure p. More precisely it has been 
shown that there exists a critical pressure pc = Pc(Qf, a, b) 



such that for < p < Pc the system admits the single 
standard phase (TOs(p), £s(p)), while a second fluid-rich 
phase (mf(p), £f(p)), appears for p > pc- The standard 
phase is similar to the unique phase described by the 
model with potential energy densitj*i^ ^I^b. In FIG.[T]the 
standard and the fluid-rich phases are depicted for p > Pc 
and for a particular choice of the physical parameters 
a, a, b. 



IV. COEXISTENCE 

The second gradient model has the same phases as the 
associated first gradient model. The main result of this 
paper is the existence of pco = ^00(0^7 a, b) > Pc, called co- 
existence pressure, such that the standard and the fluid- 
rich phases coexist at the pressure p = Pco and do not 
coexist at p > Pc and p ^ pco- The proof will be achived 
in two steps: first we shall show that there exist a unique 
value of the pressure such that the total potential en- 
ergy densities evaluated at the two phases are equal; the 
second step will consist in proving the existence of the 
connection, that is the equilibrium profile connecting the 
two phases. 



A. Coexistence pressure 

We first review some of the results in the previous pa- 
per— 1^; there we have studied the equations ^'^ = and 
^'m = looking for the minima of the function . 

We have shown that the standard phase (ms(p), £s(p)) 
is the solution of the two equations m = be and p ~ fi (£), 
for any p > 0, where fi{e) := —e — Q!6'*£'^/3. 

On the other hand the fluid-rich phase (mf(p), £f(p)) 
is the solution, with the smallest value of £, of the two 
equations m = m^{e) and p — /+(£), where 



m+(£) 



4a 
ab-^ 
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FIG. 2: Graph of the overall potential energy ^'b(p) and *I'f(p) 
for a = 0.5, & = 1, and a = 100. 



and 




/+(£): 



For e < —2/{b^/aJa) the function /+(£) is positive, di- 
verging to +00 for e — > — oo, and has a minimum at Sc 
such that /-|_(ec) = Pc] this explains why the fludized 
phase is seen only for p > p^. Moreover it has been 
proven that for any p > pc the point (mf(p), et(p)) is a 
minimum of the two variable potential energy ^'(m, e,p) 
with p fixed, while it is a saddle point for p = pc. 

We now prove the first step of the above stated co- 
existence result. For any p > pc, we let Vl's(p) := 
'i'(ms(p),es(p),p) and *f(p) := ^'(mf(p), £f(p),p) and 
prove that 



> *f(p) forp>pco 
*s(p) = *f(p) forp==pco 

< *f(p) for Pco > P > Pc 



(10) 



that is the overall potential energy density of the stan- 
dard and the fluid-rich phases are equal only at the co- 
existence pressure. This statement has been tested on 
numerical grounds, see FIG. [5] where the graphs of the 
functions 4's(p) and ^f(p) are depicted for a given set of 
physical parameters. 

In order to prove pU)) we first compute the deriva- 
tive of the two functions ^'s(j') and Vl'f(p) (with respect 
to p); by using the chain rule, and the fact that 
(ms(p), £s{p)) and (mf(p), £i{p)) are solutions of the equa- 
tions 4'm(m,e,p) = and ^'£(to, e,p) — 0, we have 
that %{p) = e^ip), mp) = Eiip), ^':{p) = e'M, and 
^'{p)^e',{p). 

Now, since £s{p) and £f{p) are negative functions of the 
pressure, we have that both ^s(p) and ^'f(p) are decreas- 
ing functions of the pressure on the interval (pc,+oo)- 
Moreover, noted that both /i and /+ are decreasing func- 
tions (of the strain) on (— (X),ec], we have that £s{p) and 
£i{p) decrease whenp increases. It then follows that £i,{p) 
and £f(p) are negative and therefore 4's(p) and ^f(p) are 
concave on the interval (pc, +oo). 



Since the two functions ^s{p) and ^{{p) are decreasing 
concave functions on the interval (pc,+oo), in order to 
prove pU)) it is sufficient to show that ^s{Pc) < ^'f(Pc) 
and 4's(p) > '^[{p) for some p sufficiently large. The 
proof of the first remark is easy: at p = Pc the two 
variable function ^'(m,e,pc) has just the two station- 
ary points^ (ms(pc),£s(Pc)) and (mf(pc), ef(Pc))- Since 
{ms{pc),£s{pc)) is a local minimum of ^'(m, e,pc), which 
tends to -|-C!0 along every direction on the plane m-e, the 
single local minimum must be the absolute minimum; 
hence, ^s(Pc) < ^f(Pc)- The second remark follows from 
the asymptotic behavior of the two functions 5's(p) and 
^'f(p); as proven in the Appendix |X1 for p — > oo we have 



0{p 



2/3-1 



(11) 



and 



^f{p) = - 



1 



1 



tP 



(1 + a62)p2 + oip) 



P 



0{p) 



(12) 



By comparing the two asymptotic formulas (jlip and ()12p 
we get immediately that for p large enough \I's(p) > 
'flip). 



B. Connection profile 

It is worth remarking that the variational problem ([3]) 
for profiles with fixed values at the end points £i and £2 
of the interval Bg, is nothing but the Hamilton princi- 
ple for a two degree of freedom mechanical system with 
Lagrangian coordinates e and m, kinetic energy T and 
potential energy U respectively given by 



Hm'f 



and 



U{m,e) = -*(m,e) 



(13) 



(14) 



and the space variable Xg interpreted as time. In other 
words the function $ defined by dS)), (O, and ([8]) is the 
Lagrangian for such a two degree of freedom equivalent 
mechanical system. 

It is important to remark that the mechanical inter- 
pretation is correct only when T is a positive definite 
quadratic form. It is easy to prove that this is the case 
provided kik^ — k^ > 0. In the limiting case fcifca — A:| = 
the form T is positive semidefinite, indeed if we substi- 
tute fci = A;|/fc3 in ([T3| the function T becomes 



1 . 



T{m',e') ^ Kim\e') = ^k^ike' + m'f (15) 



where we have set k 
ke' + m' = 0. 



fc2//c3, and is equal to zero when 
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We study now the case fcifca — fcj > and postpone the 
degenerate fci fcs — fc| = to the following section. Let us 
denote Xg by t and the derivative taken with respect to 
t by the dot. By using Q with ^ = T — U and recalling 
([TS]) , we have that the equations of motion are 



+ fci£ = — — and k^m 
oe 



dm ^ ^ 



We note that the mechanical energy of the associated 
mechanical problem E{rh,e,m,e) := T(jn,e) + U{m,e) 
is a constant of the motion. 

First note that the two points (ms(p), es(p)) and 
(mf(p), ef(p)), with p > p^, are maxima of the po- 
tential energy U{m,e) of the equivalent mechanical 
system. The problem of the existence of a connec- 
tion between the standard and the fluid-rich phase 
can be rephrased as follows: look for a solution of 
the equations (|16p . namely, a motion {mp{t),ep{t)) 
of the equivalent mechanical system, on M connect- 
ing the phase space point {mp{—oo), ep{~oo)) — 
(ms(p),£sb)) and (mp(-oo), ep(-oo)) = (0,0) to the 
phase space point (mp(-l-oo), ep(+oo)) = (mt(p), £f(p)) 
and (mp(-l-oo), ep(+oo)) = (0,0). The connection we are 
seeking for is an heteroclinic solution of the equation of 
motion tending to two fixed points in the phase space for 
t — > — cx) and t — > +oo. 

Recall that the mechanical energy E is a. constant of 
the motion and remark that at the equilibrium points it 
is equal to E{0,0,ms{p), Ss{p)) = -5'(tos(p), es(p)) and 
i?(0, 0, TOf(p), ef(p)) = — ^'(mf(p), ef(p)). From the re- 
sults in Section ITVl it follows that those two energies are 
equal only for p ~ pco- This remark yields that for any 
p > Pc and p 7^ pco the standard and the fluid-rich phases 
do not coexist. 

We are left with the case p = Pco- In principle an hete- 
roclinic solution can exist, but to prove its existence is an 
highly not trivial problem which has been solved, under 
suitable hypotheses on the potential energy, in the recent 
paper— whose main results have been summarized in the 
Appendix |B] Since in the not degenerate case the form 
T is positive definite, it is possible to find an orthogo- 
nal transformation of the coordinates in the plane m-e 
which diagonalize the form itself. Then, performing this 
transformation and subctracting to the potential energy 
U of the equivalent mechanical system the constant term 
U{msiPco),£s{Pco)) = U{m[(j)co),£i{pco)), the problem of 
finding a connection between the standard and fluid-rich 
phase is transformed in a problem in the form (|Bip with 
n = 2 and W replaced by —[U — ?7(ms(pco), es(Pco))]- 
Since this function satisfies the hypotheses of the Theo- 
rem 3.6 by Alikakos and Fusco^^ (see the Appendix [B|) 
we can then conclude that in the case p = Pco there ex- 
ists a connection between the standard and the fluid-rich 
phase and hence the two phases coexist. 




FIG. 3: Graph of the constraint curve in the plane m-e. The 
three disks represent the fluid-rich phase, the standard phase 
and the saddle (gray disk) of the potential energy Pa- 
rameters: a = 0.5, b = 1, Q = 100, fca = 1, fc = 1, and 
P = Pco = 0.24218. 



THE DEGENERATE CASE 



Consider the case k2 — iV^ife and the change of vari- 
ables 

m + ks —km + e 

and y 



where we recall k = k2/k^ = ±\/fci/fc3, which amounts 
to perform a rotation of the cartesian reference system in 
the plane m-e. Using the new variables the two functions 
T and U become respectively 

K{x, y) = T{m{x, y),e{x, y)) = ^^3(1 + k^)x^ (17) 
and 

V{x,y) = U{m{x,y),e{x,y)) (18) 
where we have used (flSll and (IT4l) with 



X — ky kx + y 

m = , and e = 



The expression ([T8|) of V is awful, but this will not be 
a problem since V is precisely the two variable functions 
— 5', which we have already deeply studied^, written via 
a rotation of the cartesian reference system. 

We apply, now, the variational principle ([3]) to the to- 
tal poroclastic potential energy density ^(i, y, x, y) := 
K{x, y) — V{x, y) and get the analogous (indeed it is a 
particularization) of the equations Q 

d^ ±d^_Q (19) 
dx dt dx dy 
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FIG. 4: Graph of the constraint curve in the plane x-y. The 
three disks represent the fluid-rich phase, the standard phase 
and the saddle (gray disk) of the potential energy Pa- 
rameters: a = 0.5, b = 1, a = 100, — 1, k = 1, and 
P = Pco = 0.24218. 



which must be solved with the boundary conditions ([5|). 
By using the definition of ^ the above equations become 



fc3(l + fc2)i 



dV 
dx 



(x, y) and 



dV_ 

dy 



{x,y)=Q (20) 



We remark that the second of the equations above is 
an algebraic equation involving the two variable x and y; 
provided it can be solved w.r.t. y, the first one becomes a 
second order ordinary differential equation in the unique 
unknown function x. More precisely, the root locus of 
dV{x,y)/dy = is made of a certain number of max- 
imal components such that each of them is the graph 
of a function x G R y{x) G M.; for each of them the 
first of the two equations (|20|) becomes a standard one di- 
mensional conservative mechanical system with potential 
energy V{x,y{x)). 



FIG. 5: Function V{x,y{x)) in the case a = 0.5, b 
100, ka = 1, = 1, and p = Pco = 0.24218. 



1, a 



Since (tos(p), £s(p)) and (TOf(p), ef(p)) satisfy the equa- 
tions 5'„i(m,e) = and 'i>e{m,e) = 0, we have that the 
two points (a;s(p), ys(p)) and {x{{p) , y{{p)) are solutions 
of the constraint equation dV{x,y)/dy ~ and hence 
they belong to the constraint curve. 

We consider, now, the case in which at p ^ pco the 
two points above fall on the same maximal component 
of the constraint equation (see FIG. [3] and S]). Using the 
conservation of the mechanical energy of the equivalent 
one dimensional conservative system allows for reducing 
the computation of the coexistence profile (heteroclinic) 
to the evaluation of a definite integral. Since the func- 
tion V has two isolated absolute maximum points which, 
by hypothesis, belong to the same maximal component of 
the constraint curve, we have that the function V{x, y{x)) 
of the real function x has two absolute isolated maxima 
in Xs{pco) and Xf(pco) (see FIG.[S]). Consider the motion 
of the equivalent one dimensional system corresponding 
to the energy level Kiax := V{xsiPco),ysiPco))- The con- 
servation of the mechanical energy implies 



A. The degenerate case: heteroclinic 

The function V is obtained by flipping the sign of the 
function 'I' and rotating the coordinate axes. This implies 
that the function V, a.t p — Pco, has the two absolute 
maximum points 



ixsip),ysip)) = ( 



nis ip) + kss (p) -krus (p) + Es (p) 



and 



{xiip),yfip)) = ( 



mf{p) + ksiip) -krriiijp) + ef{p) 



corresponding, respectively, to the standard and to the 
fluid-rich phases. 



-ksil + k')x' + Vix, yix)) = Knax 

Hence, the heteroclinic connecting the two maxima is 
given by 



t = 




fc3(l + fc2) 



Vix',yix'))_ 



■dx' 



(21) 



for any x G (min{xs(pco), a;f (p™)}, max{a;s(pco), X[{pco)}) 
and for some fixed xq in the same interval. By changing 
xq it is found a family of heteroclinic orbits which are the 
same curve up to a time translation. 

Results are depicted in the FIG. [51 The e coexistence 
profile shows a bumpiiii^ close both to the standard 
and the fluid-rich phase. This behavior is due to the 
two-dimensionality of the problem: in FIG. [7] we have 
depicted the heteroclinic on the graph of the function 
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FIG. 6: Heteroclinic (coexistence profile) in tiie case a = 0.5, b = 1, a — 100, fca = 1, fc = 1, and p — pco = 0.24218. Time 
(space in the original model) on the horizontal axis and x, y, e, and m on the vertical axis respectively in (a), (b), (c), and (d). 




0.00 



FIG. 7: Graph of the function U{m,e) and the heteroclinic 
connecting the two maxima. Parameters: a — 0.5, b = 1, 
a = 100, fca = 1, fc = 1, and p = Pco = 0.24218. 



U{m, e) = — ^'(jTi, e). From the picture it is clear that the 
optimal path climbs the two hills going around the hills 
themselves. In other words the existence of the bump in 
the connecting e-profile is due to the shape of the con- 
straint curve in the plane m-e. Since the problem has 
been reduced to the computation of the heteroclinic of 
a one dimensional conservative mechanical system in the 
X variable, it is obvious that no bump can exist in the 
x-profile. On the other hand by looking closely at the 
picture in FIG. [3] and HI it emerges that the constraint 



curve is monotonic w.r.t. y and m; this implies the mono- 
tonicity of the y and m-profiles. However, it is possible 
to find values of the parameters such that the y-profile 
presents a bump. 



In Section IV Al we have proven that in the not de- 
generate case the connecting profile does exist for any 
proper choice of the parameters. A similar result does 
not hold true in the degenerate case, indeed it is possible 
to find the connection if and only if the two maxima of 
the function U lie on the same maximal component of 
the constraint curve. We have that this is not the case 
for fc > large enough, see the dashed curve in FIG. [5] 
which is associated to the value k = 1.9. It is immediate 
to remark that the two maxima do not lie on the same 
connected component, hence in this case it is not possi- 
ble to find a connection between the fluid-rich and the 
standard phase. It is worth remarking that no evidence 
of this patologic behavior is found in the case A: < 0; 
see FIG. |9] where the constraint curve and the stationary 
point of U are depicted for a = 0.5, h = 1, a ~ 100, 
k3 ^ I, k = -0.3,-0.4,-1.0, and p ^ p^^ = 0.24218. 
This case is the most interesting one from the physical 
point of view, indeed for fc2 < the coupling between s' 
and to' is negative, hence the preferred states are such 
that the two fields s(Xs) and m{Xs) are both increasing 
or decreasing. 
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-0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 



m 

FIG. 8: The constraint curves and the points representing 
the fluid-rich phase, the standard phase and the saddle (gray 
disk) of the potential energy ^. Parameters: a — 0.5, b — 1, 
a — 100, = 1, k — 0.7, 1.9 (solid, dashed), and p — pco ~ 
0.24218. 



can be found, sec FIG.[TOl by computing the integral 



A \j2[V{x,y{x))-V{x',y{x'))] ^ > 

for any x G (min{a;, a;}, maxji, x\) with x the unique (in- 
version) point in the interval (niin{xst, a;t}, niax{xst, a^f}) 
such that V{x,y{x)) = V{x,y{x)). 

The homoclinic solution corresponding to the lowest 
maximum is often interpreted as a "critical nucleus." In 
the sense that, if a dynamic evolution would be taken into 
account, one would expect that an initial condition close 
to the critical nucleus would evolve into the standard or 
the fluid-rich phase (subcritical and supercritical behav- 
ior). This behavior depends on the size of the droplet 
in the neighborhood of t = 0. Indeed in < = the pro- 
file has the value x which, for p close to Pco, is a good 
approximation of the phase corresponding to the largest 
maximum of the function V . In other words the critical 
nucleus can be seen as a droplet of the phase correpondig 
to the smallest value of the potential energy ^ plunged 
into the other phase. 




-0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 
m 



FIG. 9: The constraint curves and the points representing 
the fluid-rich phase, the standard phase and the saddle (gray 
disk) of the potential energy . Parameters: a — 0.5, & = 1, 
a = 100, fcs = 1, fc = -0.3, -0.4, -1.0 (solid, dashed, dotted), 
and p = Pco = 0.24218. 



B. The degenerate case: homoclinic 



VI. CONCLUSIONS 

In conclusion we have studied the phase transition be- 
tween the fluid poor and rich phases in the context of 
consolidating completely fluid saturated porous media. 
A second gradient model to study the existence of such 
a transition has been proposed. Moreover, coexistence 
between the two phases at the pressure Pco, defined as 
the pressure such that the total potential energy of the 
two phases is the same, has been established. We have 
also shown that at different values of the pressure the two 
phases cannot coexist. For a particular choice of the pa- 
rameters of the model it is possible to reduce the problem 
of finding the coexistence profile to the computation of 
a definite integral. We have studied the coexistence pro- 
file for different values of the physical parameters of the 
model and shown that non-monotonic interfaces exist. 
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Consider the degenerate model and suppose that the 
pressure p is larger than pc but different from pco- Sup- 
pose that the two local minima of the potential energy 
^ lie on the same connected component of the constraint 
curve whose equation is (|20p . Consider the function 
V{x,y{x)) as in Section FV Al and note that the two lo- 
cal maxima are not equal. 

We consider the homoclinic solution corresponding to 
the lowest maximum a;. In analogy with the discussion 
of the above section, the homoclinic equilibrium profile 



Appendix A: Asymptotic behavior of potential 
energy 

In this appendix we discuss the asymptotic behavior of 
the two functions '^s{p) and ^'f(p), see Section llV Al for 
p large and, in particular, prove the equations (jlip and 
([T^ . We first note that by using ms{p) ~ bes{p), we get 
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FIG. 10: Homoclinic profile (critical nucleus) in the case a — 0.5, b = I, a = 100, ks — 1, k = 1, and p = pco — 10~^pco (left) 
and p = Pco + 10~^pco (right) with pco ~ 0.24218. Time (space in the original model) on the horizontal axis and x ((a) and 
(b)) and e ((c) and (d)) on the vertical one. 



The equation fi(e) ~ p is a cubic equation in the form 
e^ + Xe+Xp ~ 0, with A = 3/(a&'*); by Cardano's formula, 
since D := (A/3)"^ + (Ap/2)^ > 0, there exists a single real 
solution given by 



£s{p) 



1 



-Xp 



1/3 



- ^Xp - Vd' 



1/3 



1 3 



1 3 



1/3 



1 3 ^3 

3^. 



1 3 \2ll/3 

J 



By using the Taylor series (1 + x)" = X^ri^o ^"('^)-''"' 
with C„(a) = a{a — 1) ■ ■ ■ {a — n + l)/n\ being the bino- 
mial coefficient, which is convergent for — l<x<+l, it 
is not difficult to prove that es{p) = -{3/{ab^)y/^p^/^ + 
(3/a6'')2/3(l/p)i/3/3 + (9(p-5/3) for p large. By inserting 

this expression in the expansion for ^/g, we get equation 

We can perform a similar computation for ^f(p). Ac- 
counting in particular for the qualitative of e{{p), which 
tends to — oo when p is increased, we shall study the 
asymptotic behavior of m^{e) for e — )■ — oo and that of 
Efip) which is the solution of the equation /+(e) = p 
when p ^ oo. The result of this analysis will provide us 
with the asymptotic behavior of ^'f(p). First of all we 
note that 



m+(£) 



2a2 



2a 



for e — )• — oo. By using ^ we then have 



^t{p) ^pef{p) + -{l 



-ab^){e,ip)r-^ + 0[is,{p)y 



for p — > oo, where we have used that £{{p) —> — oo for 
p — > oo. The function e{{p) is implicitly defined by 
the equation f+{e) = p which is pretty complicated. 
By expanding for e — > — oo the equation becomes 
— e(l + ab^) + h{e) ~ p with h{e) a function having limit 
for e — >■ — oo. Suppose p is large enough and let Sf(j)) 
be the solution of the equation above; by the qualitative 
study we get that e[{p) — > — oo for p — > oo. It is then 
easy to show that g{p) := ef{p) — [—p/{l + a6^)] tends 
to zero for p — > oo, indeed, since Sf is the solution of the 
equation above, we have that 



dip) 



(1 + ab^)efip) + p h{ef{p)) 



1 + ab^ 



1 + ab^ 







for p ^ oo, where we have used that h{e) — !• for e ^• 
— oo and £f{p) — > — oo for p — > oo. By inserting the 
obtained expression of Sf{p) in the above expansion of 
^'f(p) we get equation (fTT 



Appendix B: General result on the existence of 
connections 

In this appendix we briefly review the main results by 
Alikakos and Fuscoi^ on the existence of connections. 
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Let W : R, with n > 1, be a C2(R") posi- 

tive function satisfying the fohowing hypotheses: (1) W 
has two distinct local minima a_ , a+ e M" such that 
W{a-) = W{a+) = 0, (2) W{u) > for any u ^ a_,a+, 

(3) liminf|„|^oo W{u) > 0, (4) there exists tq in the open 
interval (0, |a_ — a+|) such that for any ^ e R" such that 
1^1 = 1 the two maps r i-> W{a± + r^) have a strictly 
positive derivative for every r £ (0,ro). Conditions (1) 
- (3) are quite natural and physically obvious; condition 

(4) is a mild technical requirement allowing for potential 
energies with C°° contact at zeroes. 

Consider the ordinary differential equation problem 



Uxx = S/W{u) 

u(— cx)) = a_ and u(+oo) = a_|- 



(Bl) 



where w : R — > M". Solutions to the problem (|B1I) are 
known in the literature as heteroclinc motions of the 
mechanical system or connection solutions in the con- 
text of phase transitions. The Theorem 3.6 by Alikakos 
and Fuscoi^ states that, under the hypotheses discussed 
above, the problem (|B1|) admits a solution. In other 
words the theorem states the existence of a connection 



under very general and mild requirements on the poten- 
tial W. The proof of the theorem is based on a direct vari- 
ational computation. More precisely the authors prove 
the existence of a critical profile of the action functional 



Aiu) 



u{x)\^ + W{u{x)) 



dx 



on the Sobolev space Wj^'^^(R,R") of functions u : R -> 
R" such that u and its weak derivative are in L^(n,R") 
for any bounded subsets £7 C R. Such a critical profile is 
the solution of the ordinary differential equation problem 

dHU). 



Compared to the standard variational calculus, see for 
instance the paragraph 8.2 in Evan's classical booki^, the 
authors have to face the lack of compactess due to the 
infinite domain R on which the solution of the variational 
problem is defined. This problem is overcame by using 
suitable costraints that are successivley removed. It is 
also worth noting that in the Theorem 3.7 the authors 
state that the connection is a minimizer of the action 
functional A{u). 
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